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Abstract. An explicit expression for the spatial curve separating the 
region of ferroelectric order ('frozen' zone) from the disordered one 
('temperate' zone) in the six- vertex model with domain wall boundary 
i-C . conditions in its anti-ferroelectric regime is obtained. 
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1. Introduction 

The six- vertex model with domain wall boundary conditions, introduced 
in [1] and solved in [2,3], has attracted much attention over the last years, in 
particular, as an example of a system exhibiting (in an appropriate scaling 
limit) spatial phase separation phenomena. The model can be regarded as a 
nontrivial generalization of the famous domino tiling problem of the Aztec 
diamond, where the celebrated Arctic Circle phenomenon was discovered 
[4]. Among the many questions concerning this kind of effects, the shape of 
q ■ the curves separating phases is of prime interest [5-10]. 

Recently, in the series of papers [11-14], a certain progress has been 
achieved in finding an analytical expression for the curve separating the 
ferroelectrically ordered (or 'frozen') region from the disordered (or 'tem- 
perate') one. By analogy with the domino tilings, this curve is called arctic 
curve. It is to be mentioned that in the anti-ferroelectric regime one more 
phase co-exists, and, correspondingly, another separation curve emerges, 
between the regions of disorder and anti-ferroelectric order [15, 16]. Here 
we discuss only the arctic curve, for which an analytic expression will be 
provided. 

As shown in [13, 14], the arctic curve can be obtained in the parametric 
form x = x(z), y = y(z), with z £ [l,oo), as the solution of the system of 
equations 

F(z) = 0, F'(z) = 0. (1.1) 
Function F(z) = F(z; x, y) depends on x and y linearly, 

jr(z) = -g- - i-f - 2 yt " + hm (kgM*))' (1.2) 

v ' z-1 z t 2 z-2At + l N^oo N v ' 

l 
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where 



A = 



a? + b 2 - c 2 



b 

t = -. 

a 



(1.3) 



lab 



and a, b, and c are the standard weights of the six-vertex model. The func- 
tion hj\i(z) = Hm{z] A, t) appearing in the last term in (1.2) is the generating 
function for a certain boundary correlation function of the model on an N- 
by-N lattice. By construction, the solution of (1.1) provides only one of the 
four portions of the arctic curve, limited by x, y G [0, k], where k E (0, 1) is 
the value of the contact point (we assume that the lattice is scaled to the unit 
square, [0,1] x [0,1]); the remaining three portions can be easily obtained 
by exploiting the symmetries of the model (see [14] for further details). 

Derivation of formulae (1.1) and (1.2) exploits a multiple integral repre- 
sentation for certain correlation function (the so-called emptiness formation 
probability) devised to discriminate ferroelectric order and disorder, and is 
based on a strongly supported conjecture on the correspondence between 
condensation of almost all roots of the saddle-point equations and the arctic 
curve [11—14]. 

Evidently, to find an explicit formula for the arctic curve one needs 
to evaluate the last term in (1.2). In [14], this term was worked out for 
the disordered regime of the model, by extending the technique of paper 
[17]. However, that approach is unapplicable to the anti-ferroelectric regime. 
Here, in order to obtain the thermodynamic limit of function hjsr(z), we 
apply the random matrix model technique of paper [18]. This allows us to 
find an explicit form for the arctic curve of the domain-wall six- vertex model 
in the anti-ferroelectric regime as well. 



We recall that the anti- ferroelectric regime corresponds to A < —1; the 
parameter t is arbitrary nonnegative. A convenient parametrization of the 
weights in terms of the crossing parameter rj and the rapidity variable A in 
this regime is 



For convenience, we mention that this parameterization can be obtained 
from the one used in [14] (suitable for the disordered regime, see also the 
appendix) simply by changing r] ^ ^ + irj and A i — > ^ — zA. The weights in 
(2.1) are real and positive for rj and A chosen real, and obeying 77 > and 
— 77 ^ A ^ r]. 

The Izergin-Korepin formula describes the partition function of the inho- 
mogeneous six-vertex model with domain wall boundary conditions. Recall 
that the model is considered on the square lattice formed by intersection of 
N vertical and N horizontal lines. Introduce two sets of rapidity variables 
i^j}f=i an d { v k\h=li an d l e t O'jki bjk, and Cjk denote weights of the ver- 
tex being at intersection of k-th. horizontal and j-th vertical lines, obtained 



2. Izergin Korepin formula and function h^{z) 



a = sinh(r/ — A) 



b = sinh(r/ + A) 



c = sinh 2r/. 



(2.1) 
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by replacing A i-)- Xj — v k in (2-1)- The partition function of the domain- 
wall six-vertex model with these weights can be represented in terms of an 
N-by-N determinant: 



Zn(\\, • • • , Aat; ui, . . . , un) 



where 



n.f,fc=i sinh(?? - Xj + v k ) sinh(r? + Xj - v k ) 

rii^<ifc<jv sinh(A fc - Xj) $mh(vj - v k ) 
x del ,Yp{Xj-v k )\, (2.2) 



■= Sinh2f ? / 2 3) 

sinh(r/ — A) sinh(?7 + A) 

For a proof see [2,3]; another proof can be found in [11, 19]. 

A formula for the partition function of the homogeneous model can be 
found by evaluating the 'homogenous limit' in the Izergin-Korepin formula, 
namely, Ai, . . . , Ajv — > X and U\, . . . , vn — > 0. The partition function reads 



sinh(?7 — A) sinh(ry + A)] 



A'- 



z w( a, . . . , a ; o, . . . ,0) = a jj^F " [sr-'rW 

(2.4) 

The derivation simply uses Taylor expansion of entries of the rows and 
columns of the determinant [2,3]. 

Now we are ready to turn to the function h]\r(z) entering the equations 
for the arctic curve. This function is defined as the generating function 

TV 

h N (z) = Y,H^z r -\ (2.5) 

r=l 

M (r) / 

Here H y N J = H y N J {X,r]) is the one-point 'boundary' correlation function 
which gives the probability of having the sole c-weight vertex of the bound- 
ary row at the rth position. It describes a peculiarity of the configurations 
of the domain-wall six-vertex model, namely, that a boundary row always 
contains a string of a-weight vertices, next a single c-weight vertex, and 

finally a string of 6-weight vertices. 

(r) 

The correlation function admits two similar but actually different 
representations. In [19] it was computed as a determinant analogous to that 
in (2.4) but with the last column modified; this representation played an 
important role in the derivation of the arctic curve in [11-14]. 

Another representation, which is crucial for us below, is based on the 
aforementioned peculiarity of the configurations, and relates this correla- 
tion function with the Izergin-Korepin formula with one inhomogeneity [20] . 
Consider an incomplete homogeneous limit in which all rapidity variables 
but one, say Ai, tend to their homogeneous limit values, A2, ...,Ajv X 
and vi,...,un — > 0. Setting Ai = A + £, where £ is a new variable, the 
connection with the Izergin-Korepin formula, when rephrased in terms of 
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generating function (2.5), reads (see, e.g., appendix A in [14]): 



Mt(0) 



sinh(r/ — A) 



N-l 



Z JV (A + ^A,...,A;0,...,0) 
Zjv(A,...,A;0,...,0) 



sinh(r/ — A — £) 
Here the function 7(£) is given by 

sinh(r/ — A) sinh(r/ + A + £) 



7(0 := 



sinh(ry + A) sinh(r/ — A — £) 



(2.6) 



(2.7) 



Hence, the last term in (1.2) can be found by evaluating the thermodynamic 
limit of Izergin-Korepin formula with one inhomogeneity. 



3. Random matrix model formulation 

In [18] the methods of random matrix models were applied to evaluate 
the thermodynamic limit of the partition function (2.4). Here we consider 
an adaptation of the same technique to evaluate in this limit the function 
(2.6). 

The starting point of the method is to consider the Laplace transform 
for the function 93(A), 



99(A) = / e* z dm(z), (3.1) 

where dm(z) = dm(z; A, rf) is some measure, which is in the case of the anti- 
ferroelectric regime is a discrete one, and whose explicit form was provided 
in [18]. Let us define the quantity 

j / \ . x deti^fc^jy [ipjXj - v k )\ 

ijV(Ai, . . . , Ajv; V\, . . . , v N ) = — — — ryr-w r- (o.2J 

lli^<fc<iv smh(A fc - \j) smh(i/j - u k ) 

Using the properties of the determinant and the symmetry of the TV-fold 
integration measure, we can rewrite this expression in the form 

Ijv(Ai, . . . , Aat; ui, ... , vn) 



Yli^j<k^N sinh(A fc - Xj) sinh(^ - v k ) 



x m / i4& N ^ i4An dm{zi) • • • dm{ZN) > (3 - 3) 

which can be viewed as the partition function of a matrix model in a 'double 
external field' (see section 2.5.4 of [21]). 

Let us consider the homogeneous limit, namely, Ax,...,Ajv —> A and 
u±, . . . , z/jv — > of the expression above. For example, evaluating the limit 
as Xj's we get 

sinUX, - X/» ~ rriV-l., 11 W z ih ^.4) 



lim 

A 1 ,...,A iv ^Ani^<^7V sinn (Afc- Aj) nf^i! l<i<fe^iV 
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and an essentially similar expression in evaluating the limit in z/^'s. Denoting 
In '■= In(X, . . . , A; 0, . . . , 0), we readily obtain 

In = Nin^w J dmNi{Zj}) ' (3 - 5) 

with 

dm N ({ Zj }) := e x{zi+ - +ZN) ]"f (z k - Zj ) 2 dm( Zl ) . . . dm(z N ), (3.6) 

l^j<k^N 

which is in fact the one-matrix model expression of paper [ I ] for the parti- 
tion function (2.4). 

Let us now consider the case when A2, • • • , Ajv — > A, but Ai = A+£, where 
£ is arbitrary. In this case one can first expand the determinant appearing 
in the numerator along its first column, and next evaluate the limit 

deti <ife<iV [e^' Zfc ] 

hm 



x 2 ,...,\ N ^x ni<jj <fc ^Arsinh(A fc - Xj) 

(_-\\N-l \(z 1 +-+z N ) JL „ 

= — l^T- Y (-l) I_1 e S2 < TT (Zk-Zj), (3-7) 

fsinhn^-in^ -2 ?' I ^ 

(smn^j YYj=\ J- 1=1 i^j<k^N 

replacing further the sum with a contour integral (see, e.g., [22]), 

N f r& rl? 

Et-D"-'^' n <*-*>= n ^-^£ n « (z 

1=1 l4j<k^N l^j<k^N JC Llj=l\ z z 3> 

(3.8) 

where C is a simple closed counterclockwise-oriented contour enclosing all 
Zj's. As a result, for the quantity In(Q '■= In(X + £, A, . . . , A; 0, . . . , 0) we 
obtain 

* (^hW / *MW) £ u ^°_ Zj) *■ ("> 

Here and below we neglect, for simplicity, all factors not contributing to 
the logarithmic derivative of ijv(£), which is the actual quantity we are 
interested in (see equations (1.2) and (2.6)). 

When considering the large iV limit of expression of the type appearing 
in (3.5) and (3.9) in the framework of the saddle-point approximation, expe- 
rience from random matrix models suggests to rescale the 'eigenvalues' Zj, 
by a factor N, namely, Zj h-> Nzj. Indeed, to have a non-trivial saddle-point 
all terms in the logarithm of the integrand have to be of the same order. 
The solution of the saddle-point equation is encoded in the resolvent W(z), 
defined as 

1 N 1 

W(z) := lim — V (3.10) 

V ' N^oo N ^ Z- Zj 
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where now the Zj's are the solution of the set of saddle-point equations 
associated to multiple integral (3.5). 

In the case of expression (3.9), rescaling the 'eigenvalues' Zj h-» Nzj, and 
simultaneously replacing z i— > Nz, we get 

x »® a (imheF^ 7 ^ ({%}) £ nf^-^o (3 " n) 

where tojv({-Zj}) := rriN({Nzj}). The crucial point is that the set of saddle- 
point equations relative to Zj's remains the same as for the case of Ijv, and 
the corresponding solution is still encoded in W{z). We have however one 
additional saddle-point equation, relative to variable z: 

£ = WW). (3.12) 

In other words, the saddle-point for the extra integration variable is simply 
given by the functional inverse of the resolvent, evaluated at £: 

Zs. P . = W~\i). (3.13) 

Inversion relation (3.13) appears in multiple settings; see related work on 
the HCIZ integral [23] and on the connection with the R-transform in free 
probability theory [24]. 

Differentiating expression (3.11) with respect to £ we find 

i N-l J dm N ({zj}) § n T ( ^ ]x zdz 

±(logI N (Z))> = -i^_lcothC + J L (3-14) 

N N Sdfh N {{z j })§^^ ) dz 



and hence in the large N limit we obtain 



lim h\ogI N (£))' = -coth^ + z s . p .. (3.15) 

jv— >oo iv 

Recalling equation (2.6) and inversion relation (3.13) we have 

lim lAi gM7(£)) = coth(^ + A + ? ? )-cothe + ^- 1 (0. (3.16) 

iV->oo IV Oc, 

Thus the only input needed for the explicit calculation of the last term 
in (1.2) is the resolvent W{z). This has been computed explicitly both 
for disordered and anti-ferroelectric regimes of the domain-wall six-vertex 
model in [18] (see also [25,26]). In appendix we sketch the inversion of the 
resolvent for the disordered regime, showing that the result of paper [14] is 
indeed reproduced by the present method. 



4. The arctic curve in the anti-ferroelectric regime 



In the anti-ferroelectric regime the large N solution of the random ma- 
trix formulation of In gives rise to a resolvent with two cuts, separated 
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by a Douglas-Kazakov type saturated region. Recasting the result for the 
resolvent of papers [18,26] into our notations, we have 

W(z) = -^[u( V z)- Uoo ], (4.1) 

where the function u(z) is defined in terms of Jacobian elliptic functions 

sn 2 u(z) := sn2 u °°> ( 4 ' 2 ) 

with elliptic nome q = exp(— ir 2 /2r]), and K is the corresponding complete 
elliptic integral of the first kind. Quantities /3, /?', and Moo := u(z)\ z=oo are 
given by 

0= *jU- *-*ss=p u ™ = - x ' (43) 

*:=*^, (4.4) 

and we also recall that in the anti-ferroelectric regime r/ > and — r/ < A < rj, 
so that x G [0,7r/2]. 

Functional relation (4.1) can be readily inverted with the result 

z = W-HQ = /3Sn > - /3/sn 9 2 / U - - K ZM . (4.5) 
T) sn z Uoo-rj sn z [Uoo - K^/rj) 

Denoting 

a := £ (4.6) 

and substituting the values of the constants f3, j3' , and given above, and 
switching to theta functions, we find 



where 



W' 1 ^) = 2a 



a- 



d{{x)$l{x - <) - tfi(x)tf 2 (x - aO 
i?i(2x-a£)i?i(0 



%i«) i?i(a(e + A + T/))' 1 ' J 

Plugging this expression into (3.16) one obtains an explicit form of the log- 
arithmic derivative of the thermodynamic limit of function /ijv(7(£)). Note 
moreover its striking similarity with the corresponding expression for the 
disordered regime, see (A. 7). 

Now we are ready to write an explicit formula for the arctic curve in the 
anti- ferroelectric regime. As in [14], we switch from function F(z), given by 
(1.2), to function /(£), defining it as /(£) := 7' (£)-^(7(£)) • Direct calculation 
gives 

f(0 = x<p(t + \) + y<p(£-ri)-¥(Q (4.8) 
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A 
A = 




Figure 1. The arctic curve of the domain- wall six- vertex 
model for various values of A, at t = 1. 



where function <p(X) is given by (2.3) and function reads 
#(0 + A) 

= coth£ - coth(£ + A 



A)Toc^i 10gM7( ^ 



(4.9) 



tfi«) <?i (<*(£ + A + 7?)) 

Note that varying z over the interval [1, oo) corresponds to £ 6 [0, 77 — A] and 
since function 7^) is monotonously increasing over this interval we have 
7'(£) 7^ that allows us to describe the arctic curve in terms of function 
/(0 instead of F(z). 

Equations (1.1) are replaced by the equations 

no = 0, no = 0, (4.10) 

and solving this linear system for x and y we obtain: 

- vMO - - Wit) 



X 



y 



vce+AMf-^-^ce-T/Mc+A)' 
y(e+A)^(o-^(e+A)^(o 



(4.11) 



v(e + amc - ??) - ^ - r;M£ + a) ' 

We recall that here functions 93(A) and !^(£) are given by (2.3) and(4.9), 
respectively. 

Formulae (4.11) give one portion (out of four) of the arctic curve in a 
parametric form, with £ being the parameter of the curve, £ G [0,77 — A]. 
The value £ = corresponds to the contact point of the arctic curve with 
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the x-axis, and as £ increases, the whole curve is monotonously constructed, 
up to the contact point with the y-axis, at £ = r/ — A. 

The remaining three portions of the arctic curve can be easily con- 
structed using symmetry considerations, along the lines given in [14]. The 
only modification, specific of the anti-ferroelectric regime, concerns the cross- 
ing symmetry transformation (which involve exchange of the weights a -B- b), 
which for the parameterization of weights (2.1) requires changing the sign 
of A, A h-> -A. 

In figure 1 we plot (one of the four portions of) the arctic curve of the 
model for several values of A, restricting for simplicity to the case of t = 1 
(i.e., when weights obey a = b). For — 1 ^ A < 1 (disordered regime) we 
plot the expression for the arctic curve given in [14], while for A < — 1 
(anti- ferroelectric regime) we are plotting expression (4.11) above. In the 
limiting case of A — > — oo the arctic curve can be easily computed by other 
means (see, for instance, [18]). It is apparent from the figure, and can be 
easily checked analytically, that the general expression (4.11) for the arctic 
curve indeed reproduces in this case the expected result, i.e., a straight line 
segment. This is an additional indication for the validity of the 'condensation 
hypothesis' of papers [13, 14], which is at the basis of formulae (1.1) and (1.2). 
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Appendix A. Inverse resolvent for the disordered regime 

Here we sketch the derivation for the case of the disordered regime, 
showing that indeed the result of [14] for the last term in (1.2) is reproduced 
through the inversion of the resolvent. 

The disordered regime corresponds to — 1 A < 1 and the weights can 
be parameterized as follows 



a = sin(A + rj), 6 = sin(A — 77), c = sin2?7, (A.l) 
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where A and rj are real, and obey the conditions < rj ^ n/2 and rj ^ A ^ 
7r — fj. Function is defined as 

/w sin 2?? , . 

^ (A) = si^A + ^si^A-r?)- (A - 2) 

Its Laplace transform provides an integral representation which leads to the 
random matrix model formulation, now with a continuous measure with 
support on the real axis [18]. 

Repeating the procedure of Section 3 leads to the following formula for 
the last term in (1.1) in terms of the resolvent 

lim lAiog/ l7V ( 7 (^)) =co t(e + A-r ? )-cotC + W- 1 (0. (A.3) 

jV->oo iV <7£; 

The large N solution of the random matrix formulation of In in the disor- 
dered regime gives rise to a single cut resolvent (see [18,25]), given explicitly 
by 



W(z) = - + — log 
a la. 



\/ z + 2a tan n — W z tan 2 x — 2a tan x 
y^tan 2 x+1) 



(A.4) 



where now 

7r(7r - r) - A) ir 
x := —, — a : = . (A. 5) 

The functional relation W(z) = £ can be easily inverted, 

TT . i/^ „ tan 2 (x — af ) + 1 , . „. 

z = W- l {t) = -2atznx y q) . (A.6) 

tan z (x — a{J — tan z x 

An elementary calculation leads to the formula 

W~ l (£) = a cot a£ - a cot a(£ + A - rj) (A.7) 
which indeed reproduces the result obtained by other methods in [14]. 
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